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SUMMARY 


This report describes a means of extending the area of application of 
digital computers beyond the numerical data processing stage and reducing the 
need for human participation in the formulation of certain types of computer 
problems- By the use of tensor calculus and a computer language designed to 
facilitate the use of symbolic mathematical computation, a method has been 
devised whereby a digital computer can be used to do non-numeric work: that 

is, symbolic algebraic manipulation and differentiation. 

To illustrate the techniques involved, a digital computer has been used 
to derive the equations of motion of a point mass in a general orthogonal 
curvilinear coordinate system. Since this operation involves a formulation in 
terms of first- and second-order differential coefficients, it provides a good 
demonstration of a computer's capability to do non-numeric work and to assist 
in the formulation process which normally precedes the numerical data process- 
ing stage. Moreover, this particular problem serves to illustrate the advan- 
tages of the mathematical techniques employed. With the program prepared for 
this purpose the computer will derive the equations of motion in any coordi- 
nate system requested by the user. Results are presented for the following 
coordinate systems: cylindrical polar, spherical polar, oblate spheroidal 

and prolate spheroidal . 


INTRODUCTION 


Research undertaken with the object of promoting man-computer interaction 
has directed attention to the use of computers for non-numeric operations. In 
particular, the possibility of using digital computers to derive the equations 
of motion and of mathematical physics in a general curvilinear coordinate sys- 
tem has been explored. Traditionally, these functions were considered to be 
the exclusive preserve of the scientist. Nevertheless, as is shown in this 
report, digital computers can participate in the performance of such tasks. 
These computers have been used almost exclusively in the past and continue to 
be used extensively for numerical analyses of all kinds. In fact, excessive 
preoccupation with their arithmetic capabilities has tended to obscure their 
potential for non-numeric operations. Nevertheless, the extensive logic and 
storage capabilities of these computers, combined with the evolution of new 
computer languages, enable them to be used for a wide range of non-numeric 
operations. The author is aware of only two previous attempts to use comput- 
ers in this manner: (a) Reference 1 describes an interesting technique 



whereby a digital computer uses the method of Lagrange to derive equations of 
motion. The technique, as described, was not completely satisfactory in that 
part of the operation has to be performed manually, (b) Reference 2 describes 
how an IBM 709^ computer, equipped with a Formac compiler, was used to obtain 
the Christoffel symbols of the first and second kind for 12 orthogonal curvi- 
linear coordinate systems . 

If the extensive logic and storage capabilities of these computers are to 
be used to full advantage, a departure from conventional techniques of formu- 
lation may be necessary. For example, when conventional methods are used, the 
form which the equations of motion and of mathematical physics assume depends 
on the coordinate system used to describe the problem. This dependence, which 
is due to the practice of expressing vectors in terms of their physical com- 
ponents, can be removed by the simple expedient of expressing all vectors in 
terms of their tensor components. 

As a consequence of the geometrical simplification inherent in the tensor 
method, the operations involved in formulating problems in unfamiliar curvi- 
linear coordinate systems can be reduced to routine computer operations. It 
is this aspect of the tensor method which is so attractive for the types of 
computer applications contemplated. It is the purpose of this report to use 
the tensor method to show that digital computers can be used to do non-numeric 
work. With this object in mind, a computer program was written to demonstrate 
the effectiveness of the proposed technique. This program, in the Formac com- 
puter language, was used to derive the equations of motion of a point mass in 
a variety of curvilinear coordinate systems. To derive the equations of 
motion of a particle by this method, the user need only know the coordinate 
transformation equations relating the curvilinear coordinates to an orthogonal 
Cartesian triad. When this program is used and the coordinate transformation 
equations are supplied as input, the computer will derive the equations of 
motion. The equations of motion obtained will be relative to the curvilinear 
coordinate system specified by the coordinate transformation equations used as 
input. The computer presents the results in Fortran language. However, for 
the convenience of readers, the Fortran statements are translated to conven- 
tional mathematical symbolism. 


NOMENCLATURE 


A vector 

4 physical components of A 

A x (x) contravariant vector components in the x coordinate system 
Ai(x) covariant vector components in the x coordinate system 

A 1 ^(x) components of a contravariant bivector in the x coordinate system 
Aj_j(x) components of a covariant bivector in the x coordinate system 

A x (x) components of a mixed bivector in the x coordinate system 
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A i,j 

SjU) 

ai(x) 

a X (x) 

a 1 (x) 

B i (y) 

B j(y) 

B lj (y) 

B ij (y) 
B J(y) 

S j ( y ) 

b^(y) 

F ± (y) 


F jW 


g ij 

gij 

I lJ 

M 

r 


r 


T 

T 1 


covariant derivative of a contravariant vector 

covariant derivative of a covariant vector 

system of base vectors in the x coordinate system 

system of unit vectors in the directions of a^(x) 

system of base vectors reciprocal to a^(x) 

system of unit base vectors in the directions of sT^x) 

contravariant vector components in the y coordinate system 

covariant vector components in the y coordinate system 

components of a contravariant bivector in the y coordinate system 

components of a covariant bivector in the y coordinate system 

components of a mixed bivector in the y coordinate system 

system of base vectors in the y coordinate system 

system of base vectors reciprocal to bj(y) 

covariant component of gravitational force vector in the y coordi 
nate system 

covariant component of gravitational force vector in the x coordi 
nate system 



a 1 * a J 

components of the inertia tensor in the x coordinate system 

components of the inertia tensor in the y coordinate system 

mass of space vehicle 

position vector 

scalar magnitude of r 

thrust vector 

contravariant component of the thrust vector 
covariant component of the thrust vector 



U«J(x) 

V i (y) 

i 


contravariant component of velocity in the x 
contra variant component of velocity in the y 
system coordinates 


x 1 (y 1 ; y 2 >y 3 ) functional form of the transformation from the 
system to the x coordinate system 

y 1 system coordinates 

y i (x 1 ; x 2 ^x 3 ) functional form of the transformation from the 
system to the y coordinate system 


[ij >k] 


displacement along the axis of a cylinder 
Christoffel symbol of the first kind 



a 


1 

P 1 


6 1 

J 

e 


T 
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VCp 




Christoffel symbol of the second kind 

constant coefficients 

scalar magnitude of 5^ 

scalar magnitude of a 1 

Kronecker delta 

polar angle 

physical component of the thrust vector 
gravitational potential function 
gradient of gravitational potential function 
angular displacement in longitude or azimuth 


Superscripts 

a,i,j,k,2 indices of contravariance 




Subscripts 

indices of covariance 


coordinate system 
coordinate system 

y coordinate 
x coordinate 
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ANALYSIS 


Transformation Laws 

Vector components .- To facilitate the computer processing of vectors and 
bivectors, all such entities should he expressed in terms of their tensor com- 
ponents and a corresponding set of base vectors, rather than in terms of their 
physical components and a set of unit base vectors. When referred to a gen- 
eral curvilinear coordinate system, a vector A may be expressed in the fol- 
lowing alternative forms: 

A = A^a^ = A^aJ (l) 

If, in some expression, a certain index occurs twice, this means that the 
expression's to be summed with respect to that index for all admissible val- 
ues of the index, that is. 


n 





where A 1 , Aj are the tensor components of the vector A, and a^, a J are 
the corresponding systems of base vectors. In accordance with established 
convention, contravariant components will be denoted by superscripts and 
covariant components by subscripts. It is necessary to keep in mind the dis- 
tinction between contravariance and covariance because if general coordinate 
transformations are contemplated, the transformation law for the components of 
a contravariant vector denoted by superscripts differs from that for a 
covariant vector denoted by subscripts. It must be emphasized, however , that 
the covariance or contravariance of tensor components is not an intrinsic 
property of the entity under consideration. The distinction is due to the way 
in which the entity is related to its environment, the coordinate system, to 
which it is referred. For a coordinate transformation from a coordinate sys- 
tem x to a coordinate system y given by 

y 1 = y^-(x 1 ^x 2 ,x 3 ) (2) 

the transformation law for the components of a contravariant vector A 1 is 
(see appendix A and ref. 3)* 

B«j(y) = A 1 (x) (3) 

dx 1 
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where A 1 (x) are the contravariant components in the x coordinate system and 
B^(y) are the components when referred to the y coordinate system. For the 
same transformation of coordinates, other vectors, such as the gradient of a 
scalar point function, obey a different transformation law. These are the 
covariant vectors denoted by subscripts. Assuming that the coordinate trans- 
formation is reversible and one-to-one, the appropriate transformation law for 
these vector components is (see appendix A) 

Ej(y) = A t (x) (4) 

J V 

where A-^(x) are the covariant components in the x coordinate frame and 
Bj(y) are the covariant components when referred to the y coordinate frame. 
As the following argument shows , the distinction between these two transforma- 
tion laws vanishes when the transformation is orthogonal Cartesian. Let x 1 
be the components of a position vector r when referred to the x coordinate 
system which is orthogonal Cartesian. Likewise, let y J be components of the 
same vector when referred to another orthogonal Cartesian system. In this 
case the transformation of coordinates is given by 

y 1 = ot-jxj (5) 


where the cx^ are constants. The position vector r is invariant with 
respect to coordinate transformations. Hence, the square of the vector is 
also invariant. Therefore, 


ii ii i i (i k 
x J x<J = y y 1 = a^a^x u x = 6^-x x 

therefore 

a A = k 

where 6^ is the Kronecker delta, that is, (see ref. 4) 


( 6 ) 


5 


j 

k 


1 

0 


for 

for 


j = k 
j / k 


Equation (6) is the orthogonality condition which may he used to solvq equa- 
tion (5) for . If both sides of equation (5) are multiplied by 


a k. 


11,1 

a A xJ 


i i 


and 


cj J k i i 
5 k x = x = <V 


therefore , 


rJ = 


i i 


(7) 
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From equation (5 ) > it is seen that 


dxJ 


a 3 : 

J 


and from equation (7) 


8x^ i 

— - — a 

dy 1 J 


It follows from equations (8) and ( 9 ) that 

dy 1 _ dx^ 
dxJ dy 1 


(8) 


(9) 


( 10 ) 


As a consequence of equation (10 ) , the distinction between contravariant and 
covariant vectors disappears when coordinate transformations are confined to 
orthogonal Cartesian systems. This also explains why there is no preoccupa- 
tion with these vectors in the study of ordinary vector analysis. 

Ease vectors .- Subscripts assigned to a system of base vectors indicate 
that they are covariant in character and obey the covariant transformation 
law. See equation (4-) and appendix A. Therefore, if aj_(x) are a system of 
base vectors in the x coordinate system and bj(y) are the corresponding 
base vectors in the y coordinate system, then 

b; (y) ' tp 5 i (x) (11) 


In this connection it should be noted that to every system of base vector a^, 
there exists a reciprocal system of base vectors a 1 with the following 
property: 

a<j • sTf = 6^ = at • a^ (12) 

A superscript is assigned to the reciprocal base vectors to indicate their 
contravariant character, and to emphasize the fact that they obey the contra- 
variant transformation law. (See eq. ( 3 ) and appendix A.) Hence, if ai(x) 
are the reciprocal base vectors in the x coordinate system and b^(y) are the 
corresponding base vectors in the y coordinate system, then 

P( y) = ^4 a i (x) ( 13 ) 

dx 1 

In a curvilinear coordinate system the base vectors are, in general, not 
unit vectors, but are functions of the coordinates; that is, 

a-j_ = a^x 1 ,* 2 ^ 3 ) (l4) 

= aj(x 1 ,x 2 ,x 3 ) 


(15) 
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The base vectors a-j_ may be obtained as follows: _ let dr be the 
differential of a position vector r and let dx 1 be the corresponding dif- 
ferentials of the coordinates. Then by substituting dr for A , and dx 1 
for A 1 in equat ion (l) , we have 


dr = dx 

From equation (l6) the base vectors a-j_ are given by 


(16) 


dr 

dx 1 


( 17 ) 


In an orthogonal Cartesian frame of reference, the base vectors a-j_ consti- 
tute a triad of mutually orthogonal unit vectors. However, in problem formu- 
lation, it is usually convenient to use a more general curvilinear coordinate 
system. When this is done, the magnitudes of the base vectors generally 
differ from unity. 


Vector Derivatives and the Christoff el Symbols 

The scalar product of any two base vectors aq and aj may be defined as 
follows : 


a j = 


a i s ij 


(18) 


Likewise, the scalar product of the reciprocal base vectors a 1 and a^ 1 may be 
defined as 


a 1 



( 19 ) 


The symmetry of gqj and g 1 ") follows from the nature of the scalar product. 
Certain combinations of the partial derivatives of these scalar products with 
respect to the system coordinates are useful in obtaining the derivative of a 
vector, in formulating the equations of motion, or writing the equations of 
mathematical physics in a general curvilinear coordinate system. The defini- 
tions that follow are ascribed to Christoffel and are called Christoffel 
symbols (see ref. 5)* There are two of these symbols, the first of which is 
defined as 




i f d Sik d gjk 

2\ dxJ dx 1 


dx^y 


( 20 ) 


The Christoffel symbol of the second kind is 


(21) 


Derivatives of a contravariant vector .- The utility of the Christoffel 
symbols is immediately apparent when an attempt is made to find the partial 
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I 



derivative of a base vector, or its reciprocal, with respect to any system 
coordinate. Any vector A may be expressed in the forms given in equa- 
tion (l) • Furthermore, since the base vectors are, in general, functions of the 
coordinates, it follows that the derivative of A with respect to any coordi- 
nate must involve the Christoffel symbols . From equation (l) , the partial 
derivative of the contravariant form of the vector A with respect to the 
coordinate x^ is given by 

dA dA 1 _ , A i ^ a i 

- a 


Since a.^ • aj 


Likewise, 


= Sij; then 


and 


Since 


it follows that 


dap 


dx* " 

: 8x k 

L 

a Sij = 

da.- 

dx^ 

dx k 

a Sjk = 

. 

dx 1 

8x4 

8gik _ 

8a-j_ 

8x^ 

8x«5 


a i 

= — 1 

^8r N 

8x^ 

WV 


i + A v 

1 dx* 


a j + a i 


a k 


+ a j 


dr 

dx 1 


dr 


da ■ 


a& k 

8 x^ 

^ a k 

dxJ 


daj 

Sx 1 


( 22 ) 


(23) 


(2^) 


(25) 


(26) 


From equations (23) through (26) 


• a k = £x j 

8x^ 


(27) 


Therefore if equation (12) is used, the rate of change of the base vector a.^ 
with respect to xj assumes the form 

da 

chc^ 


= [ij^kja^ 


(28) 


Equation (28) gives the required rate of change of the base vector ap with 
respect to a system coordinate, in terms of the Christoffel symbol of the 
first kind and the reciprocal base vectors . A more convenient form is 
obtained if both sides of equation (28) are multiplied scalarly by the 
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reciprocal base vector a^ to yield 


^ a i —l r • • i i— k — l 

- — t • a = [ij,kja • a 


From equation (19), it is seen that 


a^ • 5^ = 


therefore 


aa i _Z 


a = [ij,k]g' 


In terms of the defining formula (2l), equation (jO) may be rewritten as 
follows: 

da i . r*l _ f U 


Therefore , 




By substitution of equation (32) in equation (22) the partial derivative of a 
vector A with respect to the system coordinate x^ is 

dA _ dA 1 _ . A i f l ] ^ 


dx^ dx' 


dA 1 _ i / l 1 


The indices i and l in the second term on the right side of equation (33) 
are dummy 1 indices, and may therefore be replaced by any other convenient 
indices, except k. To have a common base vector a^, equation (33) may "be 


rewritten as follows 


dA /dA 


,k Jk 


1 aJW = 


A J a i = A ,k a i 


Furthermore, since 


dA dx K _ dA 
dx k dt dt 


■'As already indicated, a repeated index implies summation with respect to 
that index. Since the summation index can be changed at will, it is usually 
referred to as a dummy index. Of course, the range of admissible values of 
the index must be preserved. 
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and 


dA 1 dx k _ dA^ 
^xk dt dt 


the intrinsic derivative, or the total derivative with respect to time of the 
contravariant form of the vector A, may be obtained from equation (3*0 in the 
following form: 


dA 

dt 






(35) 


where A 1 ^ is the covariant derivative of the contravariant vector A 1 with 
respect to x^. 


The natation A^j suggests that the covariant derivative of a contra- 
variant vector is not a simple covariant or contravariant vector. As the 
notation implies, A^j is a mixed tensor, with one index of contravariance 
and one index of covariance (see appendix A) . If a single- valued, reversible 
functional transformation of the form given in equation (2) is assumed, the 
transformation law for this type of entity is 


4(y> 


dy 1 dx^ a, \ 
A (x) 

dx a dy^ P 


where A^(x) are the components in the x coordinate system and B^(y) are 
the corresponding components in the y coordinate system. 


In an orthogonal Cartesian reference frame 


g 


iJ 





ij 


Therefore, since all these scalar products are constants, it follows that the 
Christoffel symbols vanish. In this case, the covariant derivative of a con- 
travariant vector reduces the sum of the partial derivatives of its physical 
components along a set of fixed axes 


dA _ dA 1 — . 
dx^ dx^ 


i = 1,2,3 


Likewise, the intrinsic derivative of a vector reduces to the ordinary time 
rates of change of the physical components along a set of fixed axes. 

For a general space of three dimensions, equation (35) assumes the form 

i - (?£ ♦ fi >* ♦ (*£ +f 2 )- + (S +f3 > (36) 


where 
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f 3 = A 


* {A} ^ ♦ {i 3 3 } *£) ♦ <{a\l ^ * U} ^ + {23} ^ 


3 l d-x 1 . J 3 l dx 2 J 3 1 dx 3 
31 f dt 132) dt 1331 dt 


The intrinsic derivative of a contravariant vector in a space of three dimen- 
sions contains 27 Christoffel symbols. However, because of the symmetry of 
the Christoffel symbols, 


and the number of independent Christoffel symbols reduces to 18. 

Derivatives of a covariant vector .- The_second alternative from equa- 
tion (1) may be used to express the vector A in terms of its covariant 
tensor components and reciprocal base vectors; that is, 

A = A i a 1 I 


In this case, the partial derivative of the vector A with respect to 


x is given by 


dA = i g-i + A . da 1 

dx^ dx^ 1 dx^ 


From equation (12) 


. a, +S 1 • = 0 

dx* 0 dx* 
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therefore 


«£ ■ a, = -S 1 • 5i 

dx* J dx k 

Substituting equation (32) in equation (43) gives 


Sai 





therefore 



( 43 ) 


(44) 


Substituting_equation (44) in equation (42) gives the partial derivative of 
the vector A with respect to in the following form 


dl_ = ^i 



(45) 


The indices i and j in the second term on the right side of equation (45) 
are dummies, and may therefore he replaced by any other indices, except k. 
In terms of the base vectors a 1 , equation (45) may be rewritten as follows 


dA 

b^ 




(46) 


where A n - v defines 

1 V 

with respect to x^. 


the covariant derivative of the covariant vector 


A i 


It may be noted that the covariant derivative of a covariant vector is 
not a vector. As the notation implies, Aj< is a doubly covariant tensor, 
that is, a tensor with two indices of covariance. If a s ingle-valued , 
reversible functional transformation, of the form given in equation (2) is 
again assumed, the transformation law for entities of this kind is 


B, 


• .(y) = — -4- 

iJ v,y . i 


■ ^4 ^4 *,„(*> 


by dy J 


up' 


where A a p(x) are the components in the x coordinate system, and Bj_j(y) are 
the corresponding components in the y coordinate system. It may be men- 
tioned in passing that moment of inertia, which is a second order tensor, has 
a transformation law of this form. (See appendix A.) 

It appears, therefore, that the operation of covariant differentiation 
of a vector or tensor increases the covariance by one index. That is, the 
xJ covariant derivative of the contravariant vector A * is A^j, which is a 
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mixed, tensor, with one index of contravariance and one index of covariance. 
The xJ covariant derivative of the covariant vector is This is 

a doubly covariant tensor or a tensor with two indices of covariance. The 
intrinsic derivative of the covariant form of the vector A is obtained from 
equation (46) in the following form 


dA = / ^i 
dt \ dt 



A j 




(47) 


For a general space of three dimensions equation (V7) assumes the form 



As in the case of the intrinsic derivative of the contravariant vector, 
the intrinsic derivative of the covariant form of the vector A is seen to 
contain 27 Christoffel symbols. However, because of the symmetry implied by 
equation (4o) , the number of independent Christoffel symbols is again reduced 
to 18. 
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Special Coordinate Systems 


The large number of terms appearing in equations (35) and (V7) is due to 
the generality of these equations which are applicable to any space of three 
dimensions. Fortunately, for the three-dimensional spaces most commonly used, 
both of these equations reduce to a more manageable form. 


For example, if base vectors of unit length are denoted by aq or a 1 , 
then in a cylindrical polar coordinate system, 


ai 

A 

= ax 

S 11 

= 1 



= x^a 2 

§22 

II 

x" 

H 

V ~rvi 

> (52) 


A 

= a3 

§33 

= 1 

j 


a 1 

= a 1 

g 11 

= 1 


a 2 

= -i- a 2 

X 1 

g 22 

1 

(xi) 2 

> (53) 

a 3 

= a 3 

§33 

= 1 

J 



As a consequence of equations (52) and (53) there are only two independent, 
nonzero Christoff el symbols in a cylindrical polar coordinate system. These 
are 



( 54 ) 


Hence, a contravariant vector referred to this coordinate system has a time 
rate of change as follows: 



(55) 


Likewise, the time rate of change of a covariant vector referred to this 
coordinate system is given by 


dA _ fdAi 
dt V dt 


A 2 




Ax 



dx 2 

dt 



CLA3 

dt 


a 3 


(56) 
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In spherical polar coordinates. 


and 


ai 

A 

= a 



Six 

a 2 

1 A 

= x-Las 


s 22 

a 3 

H 

X 

II 

sin 

x2a 3 

§33 

a 1 

= a 1 



s 11 

a 2 

II 

a 2 


g 22 

a 3 


1 

a 3 

g33 


X 1 

sin 

X 2 


( 57 ) 


= (x 1 sin x 2 )‘ 


(x 1 ) 2 


(58) 


(x 1 sin x 2 )‘ 


In this case there are six independent, nonzero Christoffel symbols, 
are 

{22} = “ xl { 3 2 3 } = " sin x2 cos x2 

{12} = {21} = ji {13} ‘ {31} = ^ 


These 


(59) 


{33} = - xl sin2 "" {23} = {3!} ’ cot * _ 

When the Christoffel symbols are substituted in equation ( 3 5).> "the time rate 
of change of a contravariant vector referred to a spherical coordinate system 
assumes the following form: 

! = (^ + a2 {*} 1 * A3 {»} *£) - 
+ + {£} ( Al ^ + a2 tO + a3 €] ^ 

+ H + {13} ( A1 ^ + A3 tO + {A} ( Aa ^ + A3 *£)] Ss (60) 

The corresponding rate of change of a covariant vector is obtained by 
substitution from equation (59) in equation (47) • In this case. 
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Alternative Derivation of the Christoffel Symbols 


In equations (20) and (21) the Christoffel symbols have been defined in 
terms of the scalar product of two of the base vectors These symbols can 
also be derived from the equations of coordinate transformation by the follow- 
ing method, which is more suitable for the applications contemplated in this 
report . 


In a rectangular Cartesian coordinate reference frame, with coordinates 
denoted by yi, an element of arc of length d's may be expressed in the fol- 
lowing form 


ds = a a dy a = a p dy P 

therefore 

ds 2 = a a • ap dy a dy P = 6^ dy a dy* 3 = dy a dy a 

Consider a curvilinear coordinate system with coordinates denoted by x 1 , and 
assume that the x and y coordinates are related by a set of transformation 
equations as follows 


y 1 = y 1 (x 1 ,x2,x 3 ) 

The element of arc d's in the x coordinate system assumes the form 


(62) 


ds = a-j_ dx 1 = a. dx^ 
J 


therefore 


ds 2 = (ai dx 1 ) • (a j dx^) = gqj dx 1 dx^ = dy a dy ( 

n a J ( 

dy dy 


.a 


a = 6y_ ^ = g _ ^i dx j 


dx 1 6x^ 


and 


r -feEh l 

W 6xi 


(63) 
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If the transformation equation (62) is reversible and one-to-one, then 

x i = x 1 (y 1 ,y 2 , y 3 ) (64) 

By substitution from equation (63) in equation (20) , the Christoffel symbol of 
the first kind assumes the following form: 

[u,k] » hi (6 5 ) 

chc^x^ dx k 

Likewise, substituting equation (63) in equation (2l) gives for the Christoffel 
symbol of the second kind 

( 66 ) 

6y a 

COMPUTER APPLICATIONS 


It is interesting to note that the preceding derivations differ from the 
conventional approach only in the method of expressing a vector in terms of 
its components and associated base vectors. Conventionally, a vector is 
expressed in terms of its physical components and a corresponding set of unit 
base vectors. The proposed method, which is the method of tensor analysis, 
uses an almost identical formulation. The important difference being that in 
this case the components are not, in general, the physical components. 

Instead, they are the tensor components which obey transformation laws corre- 
sponding to their variance. The transformation laws for contravariant and 
covariant vectors are given by equations (3) and (4), respectively. When the 
base vectors define an orthogonal Cartesian reference frame, the physical com- 
ponents and the tensor components are equal. It follows that the tensor 
method reduces to the conventional method for problems formulated in orthogonal 
Cartesian reference systems. 

As a consequence of the geometrical simplification inherent in the tensor 
method, the operations involved in obtaining derivatives and formulating the 
equations of mathematical physics in unfamiliar curvilinear coordinate systems 
are routine operations. It is this feature of the tensor method which makes 
it so attractive for computer applications . Because of their logic and stor- 
age capabilities, digital computers are well suited to such routine operations 
if they are properly programmed. 

If the Christoffel symbols are stored for all coordinate systems that are 
likely to be used, the rates of change of any vector can be obtained by a 
straightforward application of eq-uation (35) or (47). Alternatively, if the 
functional forms given by equations (62) and (64) are known, the Christoffel 
symbols may be obtained from equations (65) and (66) without using storage 
space. This technique may be illustrated by using the transformation from an 
orthogonal Cartesian reference frame to a curvilinear coordinate system in 
which are cylindrical polar coordinates. If the curvilinear coordinate 

system is cylindrical polar, the Cartesian coordinates yi are related to the 
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curvilinear coordinates by the following transformation equations (see 
sketch (a)): 



( 67 ) 


y 1 = x 1 cos x 2 
y 2 = x 1 sin x 2 
y 3 = x 3 

The inverse transformation is given by 

x 1 = Jl y 1 ) 2 + (y 2 ) 2 1 


x2 = tan-^ 


/■3 = 


= y 3 


(68) 


By substitution from equations ( 67 ) and 
(68) in equations ( 65 ) and (66) , all 
the Christoffel symbols may be obtained. 
For the special case considered there are only 2 independent, nonzero 
Christoffel symbols out of a total of l8. Of course, in practical applica- 
tions, the operation of obtaining Christoffel symbols from formulas ( 65 ) and 
(66) and the transformation equations would be performed by a computer pro- 
grammed for non-numeric operations . 


By way of illustration, equation (66) will be used to obtain the nonzero 
Christoffel symbols of the second kind. With the exception of the dummy index, 
the superscripts appearing on the right side of the equation (66) must corre- 
spond to those appearing in the Christoffel symbol . For example 

J 1 ! - d 2 y a 6xi 

^- 22 J dx 2 dx 2 8y a 

Carrying out the summation implied by the dummy index a gives 

f 1 1 c^y 1 dx 1 + d 2 y 2 dx 1 + d 2 y 3 dx 1 
dx 2 dx 2 dy 1 8 x 2 c3x 2 dy 2 dx^x 2 dy° 

Substitution for the partial differential coefficients from the functional 
relationships given by equations ( 67 ) and (68) gives 

= [-x 1 cos x 2 (cos x 2 ) - x 1 sin x 2 (sin x 2 )] = -x 1 



Likewise, 

1 = ^y 1 dx 2 + c^y 2 dx 2 + d 2 y 3 dx 2 
L 12 J dx^x 2 dy 1 dx'ddx 2 dy 2 dx^dx 2 dy 3 
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therefore 



- s - <^) 



1 _ 

xi 


The procedure for determining the Christoffel symbols for a spherical 
polar coordinate system is the same as that used for a cylindrical polar 
coordinate system. In this case, however, the terms of equation (66) have to 
be obtained from a different set of coordinate transformation equations. The 
Cartesian coordinates y 1 are related to the spherical polar coordinates xi 
by the following transformation equations. (See sketch (b).) 



(69) 


y 1 = x 1 sin x 2 cos x 3 
y 2 = xi sin x 2 sin x 3 
y 3 = x 1 cos x 2 
The inverse transformation is given by 

x 1 = ^(y 1 ) 2 + (y 2 ) 2 + (y 3 ) 2 1 


x 2 = tan 


y 1 ) 2 + (y 2 )' 
y 3 


/ (70) 


= tan-^) 


By substitution from equations (69) and (70) in equation (66) , all the 
Christoffel symbols may be obtained. For the special case being considered 
there are 6 nonzero Christoffel symbols out of a total of 18. Of course, as 
indicated previously, the operation of obtaining Christoffel symbols from 
formula (66) and the use of the transformation equations would be performed by 
a computer programmed for this kind of operation. By way of illustration 
equation (66) will again be used to obtain the nonzero Christoffel symbols. 

By substitution from equations (69) and (70) in the expanded form of equa- 

8 2 y 2 dx 1 + c^y 3 dx x \ _ 

8x 2 8x 2 (iy 2 dx 2 (3x 2 8y 3 / 


tion (66) it is found that 
f l] / 6 2 y 1 dx 1 

(22J by 1 



dx 2 

dy 1 

dx 2 

dy 1 


+ ^y 2 ^x 2 

dx 3 8x 3 dy 2 
t 6 2 y 2 dx 2 

i^X^X 2 by 2 


+ ^ 2 y 3 

c3x 3 8x 3 by 3 / 

+ JV ^ 

dx^x 2 by 3 / 


-sin x 2 cos x 2 

0 ■ i 
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/3 Way ax 3 ay d x 3 , a 2 y 3 dx 3 \ f 3 1 ■ 1 

l?- 3 j Vdx^dx 3 dy x dx^x 3 dy 2 dx^x 3 dy 3 / 131 J x 1 

J 1 ! = fJV_ + ^ + d V &cf\ = _ x i sin2 ^ 

\33j \dx 3 6x 3 dy 1 6x 3 6x 3 by 2 bx 3 bx 3 by 3 / 

J 3 \ _ / 6 2 y x dx 3 6 2 y2 dx 3 , b^y 3 dx 3 ^ / 3 \ _ A 2 

M " Wdx 3 dy 1 3x 2 3x 3 ay 2 a x 2 a x 3 a y v = W = co x 

These results are seen to agree with those obtained in equation (59) . 


The Velocity Vector 

Two methods of obtaining the Christoffel symbols have been indicated: 
one using the methods of vector calculus and the other using known differen- 
tial coefficients from the coordinate transformation equations. Since the 
latter method is more adaptable to digital logic, equations (65) and (66) 
rather than equal ions (27) and (31) are used to determine the Christoffel sym- 
bols . Given the Christoffel symbols , it is seen that there are two forms for 
the intrinsic or absolute derivative of a vector. Equation (35) gives the 
intrinsic derivative in terms of the contravariant vector components; and 
equation ( 47 ) gives the same in terms of the covariant components. Either of 
these equations may be used. However, to avoid the necessity of transforming 
covariant components into contravariant components, and vice versa, by the 
methods of appendix A, it is better to match the formula to the variance of 
the vectors. In the course of the analysis, it will become evident what the 
variance of the vectors is. For example, the variance of the differential 
elements can be determined as follows: the differential elements dy-*- in the 

y coordinate system are related to the elements dxO in the x coordinate 
system by the following equation: 

dy 1 = dxJ (71) 

oxJ 

By comparing this equation with equation (3) .> it is seen that the differential 
elements are the components of a contravariant vector. Likewise, equation ( 71 ) 
shows that the components of velocity in the y coordinate system are related 
to those in the x coordinate system by the equation 

dy 1 _ b^_ dx^ 
dt 3x j dt 

that is, 



where V 1 (y) are the velocity components in the y coordinate system, and 
U-3(x) are the velocity components in the x coordinate system. Comparison of 
equation (72) with equation (3) shows that the components of the velocity 
vector also obey the contravariant transformation law. To obtain the velocity 
vector from equation (35) , the position vector r is substituted for the 
vector A; that is. 


A = a \ ± = r (73) 

Hence, in a cylindrical polar coordinate system 

A 1 = x 1 , A 2 = 0 , A 3 = x 3 (7*0 

By substitution of these values in equation (55) , the velocity vector is 
obtained as follows 




a 3 


When the appropriate value of the Christoff el symbol is substituted from equa- 
tion (5^), the tensor components of the velocity vector are given by 


that is, 




dx 3 _ 


dt 


a 3 



(75) 


In order to reduce equation (75) to the conventional form, where the physical 
components of velocity are associated with a set of unit base vectors , equa- 
tion (Bl) may be used to express the base vectors in unitary form. In this 
form the velocity V is given by 

V = a a + (xl Sg) a 2 + ^ a 3 (76) 

dt \ dt/ dt 

If the coordinate x 1 is identified with the radial distance r, the coordi- 
nate x 2 with the polar angle 9 , and the coordinate x 3 with the axial dis- 
placement z, the equation for the velocity in a cylindrical polar coordinate 
system assumes the familiar form 


V = 


dr 

dt 


a i 


+ 




(77) 


where ax, a 2 , and a 3 are a triad of mutually orthogonal unit vectors in the 
directions of increasing r, 6 , and z, respectively. 
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In a spherical polar coordinate system, the vector r has the following 
components . 


A 1 = x 1 , A 2 = A 3 = 0 


When these values are substituted in equation (60), the velocity vector in 
this coordinate system is given by 


V = 


dx 1 — 
■dt ai 




Again, by substitution of the Christoffel symbols from equation (59) > the 
velocity vector may be expressed in terms of its tensor components and a 
corresponding set of base vectors as follows: 


that is. 


Tr dx 1 _ dx 2 _ 
V = “dt ai + “dt a2 


dx 3 _ 
+ -dt as 


v _ dx 1 — 

v - -at ai 


(78) 


From a comparison of equation (75) with (78), it is seen that when expressed 
in terms of its tensor components, the velocity vector has the same form in 
both coordinate systems. This is true, in general, since by definition 


- _ dr _ dr dx 1 
" dt ~ dx 1 dt 


(79) 


And substitution from equation (17) in equation (79) gives 


V 


dx 1 _ 
dt a i 


Of course, the physical components of velocity are different, as can be seen 
when the methods of appendix B are used to reduce the base vectors to unitary 
form. By substitution from equation (57) in equation (78), the velocity vec- 
tor may be expressed in terms of its physical components and a set of unit 
base vectors as follows: 

V - VE *1 + ( xl s “ ** vr :)® 3 (8o) 

When the coordinate x 1 is identified with the radial distance r, the coordi- 
nate x 2 with the polar angle 0, and the coordinate x 3 with the azimuth 
angle \[/, the equation for V assumes the more familiar form 

7 = I 4 l i r f)® 2 + ( rsine i ) a3 (8l) 
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where ai, a 2 , and a3 are a triad of mutually orthogonal unit vectors in the 
directions of increasing r ,8, and \|r, respectively. 


The Acceleration Vector 


If the acceleration vector were required to formulate the equations of 
motion of a particle, the velocity vector V would be substituted for the 
vector A in the equation for the intrinsic derivative. As shown in equa- 
tion (72), the velocity vector assumes the contravariant form; therefore ,_if 
the tensor components of velocity are substituted for the components of A in 
equation ( 35 )> 


A = 



(82) 


Hence, in a general curvilinear coordinate system, the acceleration vector is 
given by 


dV 

dt 



(83) 


By substitution from equation ( 82 ) in equation (83) , the acceleration vector 
may be written in the following alternative form 


dV = 
dt 



( 84 ) 


This equation gives the acceleration in any coordinate system, provided the 
Christoffel symbols are appropriate to the coordinate system chosen to 
describe the problem. 


In a three-dimensional cylindrical polar coordinate system, equation (83) 
reduces to the form given by equation ( 55 ) > when the vector V is substituted 
for the vector A. Likewise, in a three-dimensional spherical polar coordi- 
nate system, equation (83) reduces to the form given by equation (60), when 
the vector V is substituted for the vector A. If equation (60) is used to 
obtain the acceleration vector, the tensor components of velocity, rather than 
the physical components, must always be used. The tensor components of veloc- 
ity are given by equation (78) • These are 


A 1 = — 

A dt ' 


A 2 = 


dx 2 

dt 


A 3 = 


dx 3 

dt 


( 85 ) 


Substituting these values in equation (60) gives the acceleration in terms of 
spherical polar coordinates 
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dV 

dt 


+ 


f d 2 x x 
\ dt 2 

♦{£ 

1 dx 2 dx 2 
-J dt dt 

( i 2 x 2 

V dt 2 

H 

f 2 1 dx 1 dx 2 
(12J dt dt 

/d 2 x 3 

V dt 2 

H 

( 3 1 dx 1 dx 3 
[13J dt dt 



1 

1331 


dx 3 dx 3 
dt dt 


ai 


2 

33 


\ dx 3 dx 3 
>J dt dt 


a 2 


( 86 ) 


By substitution of the Christoff el symbols from equation (59) equa- 
tion (86 ) y the acceleration vector may be expressed in terms of its tensor 
components and associated base vectors as follows: 


dV 

dt 


^ - x 1 sin 2 x 2 (^ 

L dt 2 \ dty V dt 


ai 


d 2 x 2 , 2 dx 1 dx 2 . p 

+ sin xr cos 

dt 2 xi dt dt 

d 2 x 3 , 2 dx 1 dx 3 


x 2 ^' 
\ dt . 


a 2 


dt" 


. _ , 0 2 dx 2 dx 3 ^ — 

+ — r — TT" — TV + 2 cot X — TV — TV J a 3 


-i dt dt 


dt dt J 


(87) 


The corresponding physical components of the acceleration vector may be 
obtained from the base vectors expressed in terms of unit vectors in accor- 
dance with equation (57)* When appropriate substitutions are made, equa- 
tion (87) gives 


d V = 
dt 


^ - x/i 

dt 2 \ dt/ \ 


sm x 


.2 dx" 
dt 


2n 


ai 


x 


1 + 2 _ x i S in x 2 cos x 2 ! 


dt" 


dt dt 


dx£\ 


a 2 


+ ,x 2 sinx^ 


dty j 

+ 2 sin X 2 + 2x ! COB x2 Sg ag 

dt dt 


dx 2 dx 3 
dt dt 


Equations of Motion in a General Curvilinear Coordinate System 

In using tensor methods to derive equations of motion, it is again 
important to remember that the acceleration and force vectors must always be 
expressed in terms of their tensor components rather than their physical com- 
ponents. Hence, the two sides of every equation must balance with respect to 
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their covariant or contravariant properties before applying Newton's second 
law of motion. In this connection it is worth noting that, although the 
acceleration vector is expressed in contravariant form in equation ( 84 ) , the 
force vector may appear in the form of a covariant tensor. The force vector 
assumes the covariant form when it appears as the gradient of a scalar point 
function. This occurs in the equations of motion of a space vehicle which, in 
addition to the thrust force, is subject to gravitational forces. If the 
gravitational forces are expressed in the form of the gradient of a gravita- 
tional potential function, the force vector is 

F = 79 + T (88) 


where cp is the gravitational potential function, which may include the_ 
influence of oblateness and extraterrestrial gravitational forces, and T is 
the thrust vector. (See refs. 6 and 7 -) 


The tensor form of the gradient of a scalar point function assumes the 
following form 



( 89 ) 


The use of the reciprocal base vector a 1 in equation (89) is justified by 
the following considerations: the components of the gradient of the gravita- 

tional potential f met ion in the y coordinate system are related to those 
in the x coordinate system by the following equation: 

dep dep dx^ 

8y^ dx^ dy ^ 

or 

F i (y) = F ^ (x) ( 90 ) 

dy 1 J 

where 

^ ■ # 


F o< 



The transformation of the components of the gradient vector from the x coor- 
dinate system to the y coordinate system (eq. (90)) obeys the covariant 
transformation law as defined in equation ( 4 ) . 

The equation of motion of a point mass which is subject to gravitational 
and thrust forces is obtained by combining equations ( 84 ) and (88) . 
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m ("iM^ = v<p + 5 


(91) 


where M is the mass . 


It is seen that the acceleration domponents represented by the left-hand 
side of this equation are all contravariant. The thrust vector, on the other 
hand, is usually given in terms of its physical components, and as already 
indicated in equation (90) , the gravitational forces assume the form of 
covariant vectors. To have a force system compatible with the accelerations, 
it is necessary to convert all the force terms to the contravariant form. The 
potential gradient function may be converted to contravariant form with the 
aid of equation (A22) . From equations (89) and (A22) 


VCp = — — a.j = g^ — ^ a,- 

8 x j 8xJ 1 


(92) 


The thrust vector may be expressed in the following alternative forms 

T = T^ = t 1 ^ 


where T 1 are the contravariant tensor components of the thrust vector, and 
t 1 are the corresponding physical components. The physical components of the 
thrust vector are related to the contravariant tensor components by 
equation (Bl) 


T 1 = X 


7 s (ii) 


(93) 


By substitution from equations (92) and (93) in equation (91) , the equation of 
motion assumes the following form 


M 


( i 2 x 1 fi\ dx' 3 dx^_ _ /_ii 8 cp , t 1 \ - 

+ W dt dtj ai "( g 8 x j ./iTTTT r 1 


7 g (ii) 


Therefore 


M 


d a x 
dt 2 


2 \jkj dt 


dx k ' 

dt 


= g 1J 


+ 


8x^ 


■jHTi ) 


(s>4) 


When the expression for the gravitational forces is expanded in a general 
three-dimensional coordinate system, equation (9*0 becomes 


d^x 1 

M ( d - + 


dt 2 Ijk 


i \ d X j 
jkj dt 


dx*f\ _ ii i2 _3cp i3 + 

dt. /' S 8x a 3x 2 ^~ 3 


7 g (ii) 


(95) 


However, in a rectangular coordinate system 
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g 


-*- j = 0 


for 


i i J 


and 


g (ii) - i/ g(11) 

where the parentheses imply suspension of the summation convention. 

Substituting these values in equation (95) gives for orthogonal systems 

Ji I dx j dx k \ _ 1 dcp , t 1 (g 6 ) 

\jkj dt dt/ J 


M + 


dt 2 


S(ii) dx 1 ^(ii) 


Equation ( 96 ) may be rewritten as follows: 

M d 2 x i „ , , Ji i dxj dx k; \ dcp , ^ 

M + S <iD WUTj’O* -I'M T 


or in the alternative form 


M 


( d 2 x i ‘ , f dx^ dx^ dcp , i , s 


From the point of view of non-numeric computer operations, it is more expedi- 
ent to eliminate the Christoff el symbols and the metric tensors from equa- 
tion (97)- These symbols are related to the coordinate transformation 
equations by equations (65) and (63) . Substituting from equations (63) and 
(65) in equation (97) gives 


M 


dy a dy a \ 

dxTiT/ 


d 2 x i 

dt 2 


+ 


fJ: y a d£) dxJ 
^dx^dx^ dx^y ^ 


dx^ ~ 

dt 


_^p + / dy a d y~ T i 
dx 1 Jbx^ dx^ 1 ^ 


(9d) 


As indicated previously, a repeated index implies summation with respect to 
that index. An exception to this rule occurs when repeated indices are 
enclosed in parentheses. Parentheses around an index imply that the summation 
convention is to be suspended for that index. This means that for each value 
of the index i, equation (98) must be summed on a, j, and k. For example, 
when equation (98) is summed on a, it appears as follows: 

dy 1 dy 1 + dy 2 dy 2 + dy 3 dy 3 \ d 2 X 1 

dx^^ dx^^ dx^*^" dx^^ dx ^ ■*■ ^ ' dt 2 


Mi 


L NJX' 


(&Y X dy 1 
Vdxddx^ dx^ 


d 2 ^ dy 2 + d 2 y 3 dy 3 \ dxJ dx^ 
dxJdx^- dx^ dxddx^- d x-v dt dt 


= + / dy 1 dy 1 dy 2 dy 2 + dy 3 dy 3 T i 

dx-*- J dx^ 1 ) dx( i T (3x ( -*- ) dx^-*-) dx^ "*" ) dx( id 


(99) 
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The left side of this equation must also be summed on j and k. When each of 
these indices is permitted to take the values 1, 2 , in turn, equation ( 99 ) 
assumes the following form: 


Ml 


LVhc 


bflkfl + JfL-Mr + JfL Jz 

(i) dx^ 1 ) dx( 1 ^ dx^ 1 ) dxC 1 ) dx(- 


d 2 xi 


( d 2 y x dy 1 + c^y 2 dyf + d 2 / 3 dyf' \ dxi dx 1 
W^x 1 dx 1 dx^x 1 dx 1 dx^x 1 dxV ^ ^ 


x^x 2 dx 1 dx-'vix 2 dx 1 cbc^bx 2 dx 1 . 


+ 


ifzi 

.cbt-^x 3 dx 1 

d 2 y 


ifad 




d 2 v 3 


n 3N 

2Z_ 


dx-'-dx 3 dx 1 dx-'-dx 3 dx 1 / 




-s 2 2 

d JL 




n2 3 

O V 


,dx 2 dx 1 dx 1 dx 2 dx 1 dx 1 dx 2 dx 1 dx 1 , 






a a y a 


^dx 2 dx 1 

aV a , 1 


^e + 


a 2 v 3 




ax 2 ax 2 ax 1 ax 2 ax 2 ax 1 . 


i2 2 


^ax 3 ax 1 

a 2 yi ay 1 




\2 3 

a y 


ax 2 a x 3 ax 1 

dfy2 dy2 


ax 2 ax 3 ax 1 . 


.ax 3 ax 1 ax 1 ax 3 ax 1 ax 1 


^y 1 ay 1 + ay y + 


a 2 y 3 


T* 


.ax 3 ax 2 ax 1 ax 3 ax 2 ax 1 


a 2 y x ay 1 a 2 ^ ay 2 
dx 3 dx 3 ax 1 ax 3 a x 3 ax 1 ax 3 ax 3 ax 11 , 


ar 

ax 3 dx 2 ax 1 . 

ay 3 


a 2 v 3 


v~ 

dt 2 

dx 1 

dx 1 

dt 

dt 

dx 1 

dx^ 

dt 

dt 

dx 1 

dx 3 

dt 

dt 

dx 2 

dx 1 

dt 

dt 

dx 2 

dx 2 

dt 

dt 

dx 2 

, 3 

dx 

dt 

dt 

dx 3 

dx 1 

dt 

dt 

dx 3 


dt 

dt 

dx 3 

dx 3 

dt 

dt 


dy 1 dy 3 


dy2 dy2 dy^ ay 3 


= acp + . 

dx^ J bx-W a X ( iJ dx^ 1 ^ dx^ dx^ dx^ 


( 100 ) 


The form of this equation is well suited to routine non-numeric computer oper- 
ations. The large number of terms appearing in equation (100) is due to the 
generality of this equation, which is applicable to any space of three dimen- 
sions. Moreover, since this equation is applicable to any space of three 
dimensions, it may be permanently stored in the computer. Hence, to obtain 
the equations of motion in any system of coordinates, the only information 
required is the special form of equation (2) relating that system of coordi- 
nates to the orthogonal Cartesian coordinates y 1 . 
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Computer Derivations of Equations of Motion 
in Special Systems of Coordinates 


Spherical polar coordinates .- Consider a transformation of coordinates 
specifying the relation between the spherical polar coordinates x-*- and the 
orthogonal Cartesian coordinates y 1 . In this case, equation ( 2 ) becomes 
(see sketch (b)): 

y 1 = x 1 sin x 2 cos x 3 
y 2 = x 1 sin x 2 sin x 3 
y 3 = x 1 cos x 2 

These coordinate transformation equations were supplied as input to an 
IBM 709^ computer , which was programmed for non-numeric operations • 

When the computer was instructed to perform the operations involved in 
equation (lOOj , the following output was obtained in Fortran language: 

FOR I =1, 

The expression input for Y(l) is given below. 

X(l)*FMCSIN(x(2) )*FMCC0S(X(3) )$ 

FOR 1=2, 

The expression input for Y(l) is given below. 
X(1)*FMCSIN(X(2))*FMCSIN(X(3))$ 

FOR I = 3, 

The expression input for Y(l) is given below. 

X(l)*FMCC0S(X(2))$ 


FOR I = 1 

The equation for I = 1 is given below. 

M*(P(l)-R(2)**2 .0*X(l)-R(3)**2 .0*X(l)*FMCSIN(X(2) )**2 .0)$ 

=DPHl(l)+TAU(l)$ 

FOR 1=2 

The equation for I = 2 is given below. 

M*(P(2)*X(l)**2 .0+R(l) *R(2) *X(l) *2 .0-R(3)**2 .0*X(l)**2 .0*FMCSIN(X(2 ) )*FMCC0S 
( X(2 ) ) ) $=DPHI (2 ) -HTAU(2 ) *X(1 ) $ 

FOR I = 3 

The equation for I = 3 given below. 

M*(P(3)*X(l)**2.0*FMCSlN(X(2))**2.0+R(l)*R(3)*X(l)*FMCSIN(X(2) ) **2.0*2. 0+R(2) 
*R( 3 ) *X( 1 ) **2 . 0*FMCSIN( X(2 ) ) *FMCC0S( X(2 ) ) *2 .0 ) $ 

=DPHI( 3 ) +TAU( 3 ) *X(l ) *FMCSIN(X(2 ) )$ 
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Comp/load Executive 
Time on t ime , min time, min 

950.09 0.78 0.4l 

In interpreting these Fortran statements, it must he remembered that: 

p(i) ■ 1? 

In terms of conventional mathematical symbolism, these equations assume 
the following form: 


M 


r^.x/s 

o\ 2 / o\2 

£—) - x^ sin x 2 ^L.) 

dt 2 \ 

dt/ \ dt/ _ 


_ d<P + 

dx 1 


- (- 1 )" - 


dcp 


+ X^T 2 


M (^x 1 sin x 2 ^) ^ x + 2X 1 sin 2 x 2 + 2(x 1 ) 2 sin x 2 cos x 2 ^3 ^5- 

L\ / dt 2 dt dt dt dt 


= ^ + x 1 sin x 2 t 3 

dx 3 


Because of its generality, equation (100) is applicable in all coordinate 
systems. Therefore, to obtain the equations of motion in any other coordinate 
system, all that is required is to supply the computer with the appropriate 
coordinate transformation equations. 

Cylindrical polar coordinates .- As a further illustration of the procedure 
involved, consider the equations of motion in a cylindrical polar system of 
coordinates. In this case, the coordinate transformation equations are (see 
sketch (a)): 

yi = x 1 cos x 2 
y 2 = x 1 sin x 2 


y 3 = x 3 


When these coordinate transformation equations were used to evaluate the terms 
of equation (100) , the following computer output was obtained. 
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Y(l) is given below. 


FOR I = 1 

The expression input for 
X(l)*FMCC0S(X(2))$ 

FOR I = 2 

The expression input for Y(l) is given below. 
X(l)*FMCSIN(X(2))$ 

FOR I = 3 

The expression input for Y(l) is given below. 

X(3)$ 


The equation for I = 1 is given below. 

M*(P(1)-R(2)**2 -0*X(l) )$ 

=DPH I ( 1 ) +TAU ( 1 ) $ 

The equation for I = 2 is given below. 
M*(P(2)*X(l)**2.0+R(l)*R(2)*X(l)*2.0)$ 

=DPHI(2 ) +TAU(2 ) *X(1) $ 

The equation for I = 3 is given below. 

M*P(3)$ 

=DPHl( 3 ) +TAU( 3 ) $ 

Comp/load Executive 
Time on time, min t ime , min 

O37.78 1.15 0. Ill- 


Translating these equations ffom Fortran language to conventional mathematical 
symbolism yields the following: 



x i dxf. dxf\ 

dt dt J 


dx 1 


^ * 2x1 S* 


dcp 

= 


+ X^T 



Prolate spheroidal coordinates .- Another interesting system of orthogonal 
curvilinear coordinates is the prolate spheroidal coordinates. Coordinate 
surfaces are obtained by rotating a family of confocal ellipses and hyperbolas 
about their major axes. Rotating these conic sections gives rise to a system 
of prolate spheroids and hyperboloids of two sheets. A family of planes 
through the axis of rotation completes the system of orthogonal surfaces. The 
curvilinear coordinate systems generated by the curves of intersection of 
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these surfaces are useful in certain quantum mechanical problems (ref. 8). 

The transformation equations relating this system of coordinates to the 
orthogonal Cartesian system are as follows: 

yi = a sinh xi sin x 2 cos x 3 
y 2 = a sinh x 1 sin x 2 sin x 3 
y 3 = a cosh xi cos x 2 

To obtain the equations of motion relative to a prolate spheroidal system 
of coordinates, these transformation equations were substituted for equa- 
tion (2) in the computer program. Execute time was I.63 minutes. O mi tting 
the print-out in Fortran language, the equations of motion were obtained as 
follows: 

M[a 2 (sin 2 x 2 + sinh 2 x 1 ) — - X + 2a 2 sin x 2 cos x 2 
1_ dt 2 dt dt 

+ a 2 sinh x 1 cosh x 1 ^2- - a 2 sinh x 1 cosh x 1 ^2:1 ^2^. 

dt dt dt dt 


a 2 sin 2 x 2 sinh x 1 cosh xi ^2^ 

dt dt 


= a ^sin 2 x2 + sinh 2 x 1 t 1 + 


dtp 

dl 1 


f n2- r 2 

M a 2 (sin 2 x 2 + sinh 2 x 1 ) - a 2 sin x 2 cos 

L dt 2 


+ 2a 2 sinh x 1 cosh x 1 ^2_ ^2— + a 2 sin x 2 cos 

dt dt 


2 o o . , 2 1 dx 3 dx 3 

a sin x 2 cos x 2 sinh^ xi — —3 rr 

dt dt 


,0 dx 1 dx 1 

* "at It 


x 2 ^- x2 ^X 2 
dt dt 


= a (^ s in 2 x 2 + sinh 2 xi)t 2 + — ^ 

dxr 



sin 2 x 2 sinh 2 


xi 


d 2 x 3 

dt 2 


+ 2a 2 sin 2 x 2 sinh 


x 1 cosh x 1 


dx 1 dx 3 
dt dt 


+ 2a 2 


sm x 


: 2 cos x 2 sinh 2 x 1 ^2^ 


dt 


dx 3 

dt 


a sin x^ sinh x 1 t3 



Oblate spheroidal coordinates .- Confocal ellipses and hyperbolas rotated 
about their minor axes generate the oblate spheroids and hyperboloids of one 
sheet (ref. 9). These surfaces, together with a family of planes through the 
axis of rotation, constitute a family of orthogonal surfaces. The curvilinear 
coordinate systems generated by the curves of intersection of these surfaces 
are called oblate spheroidal coordinates. Oblate spheroids are sometimes 
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referred to as planetary ellipsoids, because the Earth and the planet Jupiter 
are approximately of this form. The transformation equations relating this 
system of coordinates to the orthogonal Cartesian system are as follows: 


y 1 = a cosh x 1 sin x 2 cos x 3 
y 2 = a cosh xi sin x 2 sin x 3 
y 3 = a sinh x 1 cos x 2 


These transformation equations take the place of equation (2) in the computer, 
In this case , the time required to execute the program was again 1 .63 minutes , 
Omitting the print-out in Fortran language, the equations of motion relative 
to a system of oblate spheroidal coordinates were obtained in the following 
form: 

M a 2 (sinh 2 xi + cos 2 x 2 ) + a 2 (sinh x 1 cosh x 1 ) ^2- 

dt 2 dt dt 


2a 2 cos x 2 sin x 2 - a 2 sinh x 1 cosh x 1 5^- 

dt dt dt dt 


a 2 cosh x 1 sinh x 1 sin 2 x 2 = a(\/sinh 2 x 1 + cos 2 x 2 )t 1 + -flS? 

dt dt J j x i 


M a 2 (sinh 2 x 1 + cos 2 x 2 ) ^ 2 - 2 + a 2 sin x 2 cos x 2 

dt 2 dt dt 


+ 2a 2 sinh x 1 cosh x 1 ■ - a 2 sin x 2 cos x 2 

dt dt dt dt 


- a 2 cosh 2 x 1 sin x2 cos x 2 = a-Vsinh 2 x 1 + cos 2 x 2 t 2 + — ^ 

dt dt J J x 2 


M a 2 cosh 2 x 1 sin 2 x 2 ^ — + 2a 2 sinh x 1 cosh x 1 sin 2 x 2 ^2- 

dt 2 dt dt 


+ 2a 2 cosh 2 x 1 sin x 2 cos x2 = a cosh x 1 sin x^t 3 + — ^ 

dt dt 3x3 


CONCLUSIONS 


The extensive logic and storage capabilities of digital computers, com- 
bined with the new computer languages, enable them to be used for a wide range 
of non-numeric operations. Research indicates that these computers can be 
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used more effectively for this purpose if all vector quantities are expressed 
in terms of their tensor components, rather than in terms of their physical 
components. Because of the geometrical simplification inherent in the tensor 
method, the formulation of problems in curvilinear coordinate systems can be 
reduced to routine computer operations. The results obtained suggest that the 
exploitation and extension of these techniques should lead to a substantial 
reduction in the man hours required to formulate and process engineering and 
scientific problems. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Dec. 8, 1966 
125 - 19 - 01-37 
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APPENDIX A 


TRANSFORMATION FORMULAS FOR VECTORS AND BIVECTORS 


Base Vectors 

The transformation laws and, hence, the covariant and contravariant 
character of the base vectors and their reciprocals may be obtained as fol- 
lows: Let the differential of a position vector be denoted by_ dr. Then if 

ai(x) are the base vectors in the x coordinate system, and bj(y) the base 
vectors in the y coordinate system, the differential dr may be expressed 
in the following alternative forms: 

dr = aiCxjdx 1 = b 1 (y)dyj = b^y) dx 1 (Al) 

J J dx 1 

Therefore 

=!<*) * b,(y) (A2) 

dx 1 J 

Likewise , 

a-r(x) — dy J = b i (y)dyj 

1 dyJ J 

therefore 

tj(y) = ^4 ai(x) (Aj) 

dyd 

It is seen from equations (A2) and (A 3 ) that the base vectors aq and bj obey 
the covariant transformation law; consequently, the use of subscripts is 
justified. 


Reciprocal Base Vectors 
To each system of base vectors ad there exists a reciprocal system of 


vectors with the following property 


a J = &^ = aj • a-j_ 


(AM 


where S!r is the Kronecker delta; that is, 


= 1 


for 

for 


j = i 
j t i 


Scalar multiplication of each side of equation (A2) by b J (y) gives on using 


b'Ky) • a ± (x) = 

dx. 


(A5) 


Similarly , from equations (A3) and (a 4) it is seen that 

ax 1 

dyJ 


a X (x) • b .(y) = 


(A6) 


Equation (Al) referred to the reciprocal system of base vectors assumes the 
form 


therefore 


and 


therefore 


dr = 

a 1 (x)dxi = 

= b J (y)dyj 

dy j 

= ^(y) • 

a x (x)dx^ 

dyj 

= Mdx ± 

OyJ 1 


dx-j_ 

= a ± (x) • 

^(y)dyj 

dXj^ 

= dy- 

Ox 1 J 



From equations (A7) and (A8) 


Sx 1 


therefore 


i 1 (x)dx 1 = b J (y) -=-r dxq 

oy J 


a 1 (x) = b ^(y) 


dyJ 


Likewise, from equations (A7) and (A 9 ) 

bJ(y) = a 1 (x) 

dx 1 


(AT) 


(a8) 


(A9) 


(A10) 


(All) 


From equations (AlO) and (All), it is seen that the reciprocal base vectors 
obey the contravariant law of transformation; therefore, the superscript 
notation is justified. 


Vector Transformations 

Equations (AlO ) and (All) may be used to obtain the transformation law 
for a vector A, where 
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(A12) 


A = A i a i = A^a^ 

If A = A i (x)a^(x) when the vector A is referred to the x coordinate sys- 
tem, and if A = B J (y)bj(y) when referred to the y coordinate system, 
invariance of A requires that 

B*Hy)t>j(y) = A 1 (x)a i (x) (A13) 

From equations (A2) and (Alj), the appropriate transformation law is obtained 
as follows: 

B j(y) = A 1 (x) (Al4) 

dx 1 

Equation (Al4) is the contravariant transformation law for the components of 
the vector A. When A is referred to the x coordinate system with base 
vectors aq(x) , which obey the covariant transformation law , the components 
A x (x) obey the contravariant transformation law; hence , the use of super- . 
scripts is justified. If A is referred to the reciprocal base system a , 
then from equation (A12) 


A a Aqa 1 

On a transformation of coordinates__from the x coordinate system to the y 
coordinate system, invariance of A requires that 

A i (x)a 1 (x) = Bj (y)b^ (y) (A15) 

From equations (A10) and (A15) > the appropriate transformation law is obtained 
as follows: 

Bj(y) = A i^ x ^ (A 16 ) 

oy J 

It is seen that when a vector A is referred to a coordinate system with 
reciprocal base vectors, which obey the contravariant law, the corresponding 
components of A obey the covariant law, and the use of subscripts is 
therefore just if ied . 


Raising and Lowering of Indices 

Lowering indices ■ - The vector A may be expressed in the alternative 
forms given in equation (A12) . Scalar multiplication of each side of equa- 
tion (A12) by aj gives 

(a: • a.)A^ = A-(aJ * a^) (A17) 

± J J J 


3*3 



By substitution from equations (l8) and (A4) in equation (A17 ) > the following 
result is obtained 

gijA 1 = Aj (A18) 

Again by substitution for Aj from equation (Al8) in equation (A12) 

= §ij^ (A!9) 

Equation (Al8) gives the transformation from the contravariant components to 
the covariant components of a vector. The corresponding transformation of 
base vectors is given by equation (A19) • These operations are usually 
referred to as lowering the index. 

Raising indices .- Scalar mult iplicat ion of each side of equation (A12) by 
a 1 gives 


A^a^ • a^) = Aj(a^ • a*) 

Substitution from equations (19) and (a 4) in equation (A20) gives 

A 1 = g^Aj 


(A20) 


(A21) 


When this expression for A 1 is substituted in the left-hand side of equa- 
tion (A12), the following result is obtained 

g-*-jAiai = A.aJ 
o J 

therefore 

(A22) 

Equation (A2l) enables the contravariant components of a vector to be 
expressed in terms of its covariant components* Equation (A22) gives the 
corresponding transformation of base vectors. These operations are usually 
referred to as raising the index. 


Bivector Transformations 

A second-order tensor is characterized by having two indices. Both 
indices can be superscripts, in which case the tensor is doubly contravariant. 
Tensors of this kind are sometimes referred to as the contravariant components 
of a bivector (ref. 10) . When both indices are subscripts, the tensors are 
doubly covariant, or the components of a covariant bivector. It sometimes 
happens that one of the indices is a superscript and the other one a subscript. 
Entities of this kind are called mixed tensors or the components of a mixed 
bivector . 
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Contra variant bivectors .- As in the case of vectors or first-order 
tensors, bivectors are entities whose properties are independent of the ref- 
erence frames used to describe them* Equations (A13) and (A15) are mathemati- 
cal expressions of this statement, insofar as it applies to vectors. As might 
be expected, the invariance of a bivector in going from a coordinate system x 
with base vectors aj(x), to a coordinate system y with base vectors bj_(y), 

involves the equality of two dyadics. The coefficients of the individual 
dyads in these dyadics are the components of the bivectors. If in the x 
coordinate system with base vectors aj the bivector is given by 

A af 3 (x)a a (x)ap(x) 


and if in the y coordinate system with base vectors b^ this bivector 
assumes the form 


B 1 j(y)'b i (y)bj(y) 


invariance requires that 

B 1 ^(y)b i (y)bj(y) = A aP (x)a a (x)ap(x) (A23) 

By substitution from equation (A2) in equation (A23) 

B lJ (y)b i (yJbj(y) = A ap (x) b^y) bj(y) (A24) 

Therefore, by equating coefficients of like dyads in equation (A24) , the 
required transformation law is obtained as follows: 

B iJ (y) =idlfA a P(x) (A25) 

dx a 

This is the transformation law for the components of a contravariant bivector. 

Covariant bivectors .- Since covariant bivectors are characterized by two 
subscripts, it follows that the formulation of the dyadics will be in terms 
of the reciprocal base vectors. That is, if A xx p(x) are the components of the 
covariant bivector in the x coordinate system, and Bij(y) are the corre- 
sponding components in the y coordinate system, invariance of the bivectors 
requires that 

B ij(y)b 1 (y)b^(y) = A CL p(x)a x (x)a ( 3 (x) (A26) 

Substitution from equation (A 10 ) in equation (A26) gives 

B i i (y)B 1 (y)bj(y) = A o(x) b 1 (y) b*Hy) (A 27 ) 

J P dy 1 dyJ 
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therefore 


= TT ~T Aap(x) (A28) 

Equation (A28) is the transformation law for the components of a covariant 
bivector . 

Mixed bivectors .- A mixed bivector has one index of covariance and one 
index of contra variance . In this case, the bivectors consist of base vectors 
and reciprocal base vectors . The invariance requirements may be stated as 
follows: 

B^(y)b i (y)b J (y) = Ap(x)a a (x)a P (x) (A29) 

Substitution from equations (A2) and (A10) in equation (A29) gives 

B^(y)bi(y)bj(y) = A®(x) b-^y) ^ t-i(y) 

J H ox oy d 

Therefore, 

B* 

J 

The components of mixed bivectors transform, according to equation (A30) . 


Moments and Products of Inertia as the Components of a Bivector 

Moments and products of inertia provide a good illustration of the com- 
ponents of a bivector. Moreover, since the moments and products of inertia 
are the components of a Cartesian bivector, either equation (A25) or (A28) may 
be used to transform the components . The equality of covariant and contra- 
variant components of Cartesian bivectors follows from equation (10) . 

Consider the rigid -body shown in sketch (c) . Let m be the mass of a 
particle of the body at the point P, and r its position vector relative to 

the point 0, the center of mass of 
the body. The moment of inertia I 
of the body about an axis through 0, 
in the direction of the unit vector p 
is given by 

I = Dn(rxp) • (rxp) = p • ® • p 

where 0 is the inertia dyadic which 
is defined as follows: 

Sketch (c) $ = Em(r 2 I - rr) 



(y) = V; ^ A a (x 
8x a dyJ 


(A30) 
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and I is the identity dyadic; therefore. 


® ^X 1 X ia i^'l ■*-X 1 X 2a i a 2 lx 3 ' 3 " l a 3 ^X 2 X l a 2 a 1 + ^x 2 X^ a 2 a 2 2 x 2 x 3 a a a 3 


“ I X 3 X ia 3 a l ” I x3x2 a 3 a 2 + 


I x 3x3 a 3 a 3 ) 


The components I x i x j are the moments and products of inertia. 


In order to take advantage of the tensor notation, the inertia tensor will 
he denoted by when it is referred to the x coordinate system, and 

by T i 0b i bj when referred to the y coordinate system. See sketch (d) . The 

components of the inertia tensor are 
related to the components of the iner- 
tia dyadic ® as follows: 



l(JJ) ■ l M)M) 


i ly - -i i 


X X 1 


^ j 


(A31) 


Consider a transformation of coor- 
dinates in which the transformation 
equations are: 

y 1 = x 1 cos 0 - x 3 sin 0 


y 2 = x 2 

y 3 = x 1 sin 0 + x 3 cos 0 




(A32) 


Sketch (d) 


These equations represent a rotation 
about the a. 2 axis, the angle of rota- 
tion being 0 . 


Equation (A25) may be used to find the moments and products of inertia 
relative to the y coordinate frame. These are: 


jij = fV; ^y£ jO-3 
dx a dx * 3 


(A33) 


When equation (A32) was supplied as input to a digital computer, which was 
programmed to determine the inertia components T 1 ^ according to the trans- 
formation law (A33) , the following output was obtained in Fortran language: 


b-2 



^3 


THE EXPRESSION OUTPUT FOR IBARIJ IS GIVEN BELOW 



In terms of conventional mathematical symbolism, these equations assume the 
following form: 

V y i = ° 20 W + C0S0 W + cesei X 3 xl + s 2 ei x3x 3 

I yly2 = _Cel xlx2 + S0I x3 X 2 = ^yl 

Iyly3 = C0S0I X 1 X 1 " C20I X 1 X 3 + ^^X 1 ’ 8000 = 

Vy 2 = J x 2 x 2 

I y2y3 = “ 00 ^x 3 = V^ 2 

Iy3y3 = S'" 0 W - - C9Se W + C20 I x3x3 


where 



APPENDIX B 


PHYSICAL COMPONENTS 


The transformation from covariant to contravariant components and vice 
versa was discussed in appendix A. This appendix is concerned with the trans- 
formation from covariant and contravariant components to physical components 
and vice versa (ref. ll) . 

It frequently happens that an analysis can be performed and the results 
obtained, without reference to physical components. However, sometimes a 
quantity, such as a force, is known only in terms of its physical components. 
In this case, the transformation from physical components to tensor components 
must be determined. The appropriate transformations may be obtained as fol- 
lows: Prom equation (l8) 


a i ' a «j ~ Sij 

therefore 


= S (ii) 


Let 


where as a scalar magnitude and is a unit vector. With this notation 

a i a i ( a i^ = ®(ii) 

Therefore , 

a i 

that is, 



(Bl) 


where the parentheses imply suspension of the summation convention. Hence, if 
A-*- are the contravariant tensor components of a vector A, and if 1 are 
the corresponding physical components, then 





^ ~ \/ s (ii) Al 


therefore 


(B2) 

1^5 


Likewise , from equation (19) , 


Let 


a^- • = g^ 


S 1 = ^a 1 


where (3 1 is a scalar magnitude and aA is a unit reciprocal base vector. 
Therefore , 

. 2 


a 1 • ai = g( ii} = (P 1 )' 


therefore 


3 1 


therefore 


— i / (ii) s i 
a = yg a x 


(B3) 


Hence , if Aj_ are the covariant components of a vector A, and if are 

the corresponding physical components 


A i al = (/ g (il) A i) al =«^i al 


therefore 


A i 


Moreover, if the coordinate system is orthogonal, the physical components 
can be expressed in the following alternative forms 


= Jgfc ~ A ± = — Aj_ (B4) 

■fil iiy 

Equation (Al8) may be used to show that Aj_ = A^ in orthogonal coordinate 
systems. From equations (b 4) and (Al8) 


cA ± - 


A i = 


7s(iry v s (ii) L 


sun 41 


-Tiul) a1 -* 1 


That is 




That in nonorthogonal coordinate systems may be seen as follows: 
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Therefore, in this case. 


The fact that f 
of the relation 


= A i = Sij A ^ 

^ 7 S (ii) Al = ^ 

in nonorthogonal coordinate systems is a consequence 

A i = §i,i Aj 


4t 



APPENDIX C 


TRANSFORMATION FORMULAS FOR COVARIANT DERIVATIVES 


The transformation from covariant tensor components to contravariant 
tensor components and vice versa, was discussed in appendix A- The relation- 
ship between physical components and tensor components was derived in appen- 
dix B. In this appendix, it will be shown that the covariant derivative of a 
covariant vector may be obtained from the covariant derivative of a contra- 
variant vector and vice versa . The method is analogous to the procedure of 
raising and lowering of indices, which was discussed in appendix A. 


From equations (3*0 and (46) 


( B + {j*} aJ ) Sl ’ (B 


A} a j) sl 


(Cl) 


that is 


From equation (A19) 


A i,k ai = A j,k aJ = A Vi 


'• = a 


Sij a ° = a i 


Substitution from (C3) in (C2) gives 


(C2) 


(C3) 


therefore 


that is 


A 4(«io aj) - A j,k aj 


>ij A ,k A j,k 




dx^ 




Similarly, from equation (A22) 

g 3 ^ = SO 

Substitution from equation (c6) in equation (02) gives 

A o,k (glfe i ) = A )k s i 


'a 


(c4) 


(C5) 


(c6) 


46 




49 



REFERENCES 


1 . J irauch , D. H.j and Westerwick, R. A.: A Use of Digital Computers for 

Algebraic Computation. RTD-TDR- 6 3-4200, March 1964. 

2 . Walton, John J.: Formula Manipulation Code for the Calculation of 

Christoffel Symbols. UCRL-l4l77, 5 May 1965 . 

3* Sokolnikof f , I. S.: Tensor Analysis: Theory and Applications. John 

Wiley and Sons, Inc., 1951, p. 62. 

4. Spain, B.: Tensor Calculus. Third ed., Oliver and Boyd, Edinburgh, i 960 , 

p. 4. 

5* McConnell, A. J.: Application of Tensor Analysis. Dover Pub., Inc., 

1957, P- 141. 

6 . Howard, James C.: Theoretical Study of Orbital Element Control With 

Potential Applications to Manned Space Missions. NASA TN D-2237; 1964. 

7- Howard, James C.: Application of Intrinsic Differentiation to Orbital 

Problems Involving Curvilinear Coordinate Systems. NASA TN D- 286 l, 

1965 . 

8 . Margenan, H.; and Murphy, G. M.: The Mathematics of Physics and Chemis- 

try. D. Van Nostrand Co., Inc., pp. 5-7- 

9* Korn, G. A.; and Korn, T. M-: Mathematical Handbook for Scientists and 

Engineers. McGraw-Hill Book Co., Inc., 1961 , p. 177 . 

10. Wills, A. P.: Vector Analysis With an Introduction to Tensor Analysis. 

Dover Pub., Inc., 1958 , p. 221. 

11. Dorgelo, H. B.; and Schouten, J. A.: On Unities and Dimensions. 

Koninklijke Nederlandsche Van Wetenschappen . Proc . Section of Sciences, 
Vol. XLIX (Nos. 1-5). North Holland Pub. Co., Amsterdam, 1946, p. 124. 


50 


NASA -Langley, 1967 8 A-2519 



"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

—National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 
important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless of 
importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distribu- 
tion because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Scientific and technical information generated 

under a NASA contract or grant and considered an important contribution to 
existing knowledge. 

TECHNICAL TRANSLATIONS: Information published in a foreign 

language considered to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information derived from or of value to NASA 

activities. Publications include conference proceedings, monographs, data 
compilations, handbooks, sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION PUBLICATIONS: Information on tech- 

nology used by NASA that may be of particular interest in commercial and other 
non-aerospace applications. Publications include Tech Briefs, Technology 
Utilization Reports and Notes, and Technology Surveys. 


Details on the availability of these publications may 


be obtained from: 


SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



